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Abstract A stationary wave pattern occurring in a flow of a two-component Bose-Einstein conden- 
sate past an obstacle is studied. We consider the general case of unequal velocities of two superfluid 
components. The Landau criterium applied to the two-component system determines a certain region 
in the velocity space in which superfluidity may take place. Stationary waves arise out of this region, 
but under the additional condition that the relative velocity of the components does not exceed some 
critical value. Under increase of the relative velocity the spectrum of the excitations becomes complex 
valued and the stationary wave pattern is broken. In case of equal velocities two sets of stationary 
waves that correspond to the lower and the upper Bogolyubov mode can arise. If one component flows 
and the other is at rest only one set of waves may emerge. Two or even three interfere sets of waves 
may arise if the velocities approximately of equal value and the angle between the velocities is close 
to 7r/2. In two latter cases the stationary waves correspond to the lower mode and the densities of the 
components oscillate out-of-phase. The ratio of amplitudes of the components in the stationary waves 
is computed. This quantity depends on the relative velocity, is different for different sets of waves, and 
varies along the crests of the waves. For the cases where two or three waves interfere the density images 
are obtained. 

Keywords stationary waves • two-component Bose-Einstein condensate • supersonic flow • Landau 
criterium 

PACS 03.75.Kk ■ 03.75.Mn ■ 67.85.De ■ 67.85.Fg 
1 Introduction 

A unique feature of two-component Bose-Einstein condensates (BEC) is the possibility for two su- 
perfluids to flow with different velocities. The properties of such systems can be described by the 
three- velocity hydrodynamics (one normal and two superfluid velocities) . This feature was already no- 
ticed by Khalatnikov [1]. The modern three- velocity superfluid hydrodynamic theory was formulated 
in the paper by Andreev and Bashkin [3j . As was shown in [2] the specifics of the three- velocity super- 
fluid hydrodynamics is the presence of a non-dissipative drag between the components. The drag effect 
emerges at nonzero relative velocity of the components and consists in a dependence of the superfluid 
current of one component on the gradient of the phase of the order parameter of the other component. 
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The microscopic theory of the non-dissipative drag effect was developed in (3j|4]. An important ques- 
tion that arises in the three-velocity superfluid hydrodynamics is the question on critical velocities. 
At equal velocities of the components the answer to this question can be obtained from the Galilean 
transformation. It yields that under neglecting of vortex excitations the critical velocity coincides with 
the minimal phase velocity of the lower Bogolyubov mode. In a wider sense one can introduce two 
critical velocities c_ and c+, one is for the lower mode and the other is for the upper mode [5,6J. 

In case of unequal velocities of the components the situation becomes more complicated. The 
question was considered in Refs. [71(5] • The authors of [7] have obtained the dispersion equation for 
the spectrum of elementary excitations in the presence of superfluid flows. However the analysis of 
critical velocities in [Jj was based on the implication that the Landau criterium can be formulated as 
the condition on the relative velocity of the superfluid components. Such an implication can be put 
in question since there are two independent relative velocities in the three-velocity theory. In [8] the 
Landau criterium was formulated as the condition of positiveness of energies of elementary excitations 
in a reference frame connected with a normal component. It yields a joint condition on absolute 
values of the superfluid velocities of the components and on the angle between their directions. If one 
component is at rest the superfluidity condition [8] is reduced to [7j. 

The analysis carried out in [H] shows that the Landau criterium may be fulfilled if one or even 
both superfluid velocities exceed the velocity of the lower mode c_(in the latter case the velocities 
should have different directions). In view of unusual behaviour of critical velocities in such systems 
it is interesting to consider how this behavior can reveal itself in experiments. Two-component BECs 
have been realized experimentally in ultracold alkali metals gases confined in magnetic and magneto- 
optical traps. Two components may correspond to different hyperfine Zeeman states of the same 
isotope 9, 10J, or to different isotopes pT l lT2 l [T3]. One of the methods to determine critical velocities 
for trapped ultracold gases consists in the observation of density excitations induced by some object 
moving through the condensate [14| (usually a laser beam is used as such an object). A motion of an 
object in a two-component gas with nonzero relative velocity of the superfluid components corresponds 
to the general case of the three- velocity hydrodynamics. 

The Bogolyubov spectrum has a dispersion. Therefore, a motion of an object through a superfluid 
system (or a superfluid flow past an obstacle) can lead to an occurrence of stationary waves (the 
waves whose crests remain at rest relative to the obstacle). Such an effect called "ship waves" is well- 
known |15j . It was considered by Kelvin for the waves generated on a water surface by a ship moving 
in a deep water. Stationary waves in a one-component quasi-two-dimensional BEC were studied in 
|16j . It was shown that in a superfluid that flows past a point obstacle (the obstacle size is less than 
the healing length) the stationary wave pattern is similar to one for stationary capillary waves. The 
effect takes place if the superfluid velocity s exceeds the minimal phase velocity of the Bogolyubov 
mode Co- The stationary waves arise outside the Mach cone bounded by arms directed at the angles 
■d = ± arcsin(co / s) relative to the flow. 

In recent papers the solitons [61H7] and stationary waves [TjJ induced by an obstacle in a two- 
component superfluid system were studied (the paper [T7J was published as the electronic preprint 
when the present study was almost completed). But the authors of [61IT7] considered only the case of 
equal superfluid velocities (a relative velocity of the components is equal to zero). 

In the present paper we put emphasis on the general case of nonzero relative velocity. In Sec. [2] we 
obtain the equation for the spectrum and the eigenvectors of collective excitations and define two more 
critical velocities (in addition to c_ and c+). One of them is the maximum critical velocity c m for a given 
component (c TO > c_). The superflow at the velocity s-tc ra can be reached if the other component is 
at rest. The other is the relative critical velocity c scp (c scp > c rn ). If the relative velocity exceeds c sep , 
the frequency of the lower mode becomes complex valued. The latter signals for an instability of the 
two-component system with respect to a spatial separation of the components. In Sec|3]the equation 
that describes the stationary wave pattern is obtained. It is shown that the stationary waves emerge 
if the Landau criterium of superfluidity is violated (the energies of the excitations with certain wave 
vectors becomes negative), but the system remains stable with respect to the spatial separation (the 
critical relative velocity is not achieved). A number of stationary wave patterns are presented. It is 
shown that depending on the velocities of the components several qualitatively different situations are 
possible. If the velocities are the same in modulus and in direction, one set of stationary waves appears 
at s > c_, and another set adds at s > c+. The phase separation does not occur. If only one component 
moves with the velocity s, one set of stationary waves is formed at s > c m but if s reaches c SC p the 
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phase separation occurs. If the velocities are equal in modulus and the angle between their directions is 
7r/2 or close to 7r/2 two or three sets of stationary waves occur at s\ = s 2 > c m . In Sec.|4]we investigate 
the structure of the stationary waves. The densities of the components always oscillate out-of-phase 
in the waves that correspond to the lower mode . The ratio of the amplitudes of the oscillations of 
the components depends on the relative velocity and varies along the crests of the waves. For complex 
density patterns where two or three waves interfere the density plots are presented. It is established 
that the stationary waves are visible in total density images as well as relative density images, but in 
most cases relative density images are more contrast. The only exception is the stationary waves that 
are exited at s > c + and correspond to the upper Bogolyubov mode. 



2 The spectrum of collective modes 

To analyze the stationary waves in a two-component supcrfluid system one should obtain the collective 
modes spectrum in a moving condensate. It can be found from the matrix version of the Gross-Pitaevskii 
equation 

dip 

^—^r = ~o — A^i +7i|V'i| 2 V'i +7i2|^2|Vi, 
ot 2mi 

dip Tt^ 

ifr-^- = -7, +72|V ; 2| 2 V'2 +7l2|V ; l| 2 V ; 2, (1) 

Ot 2ni2 

where ipi are the wave functions of the components, rrn are the masses of the particles, 

4:Trh 2 ai 2-Kh 2 (mi + m, 2 )ai2 

7i = . 7i2 = (2) 

nrii m\mi 

are the interaction constants (a* and a 12 are scattering lengthes). 

Here we restrict the consideration by the most convenient for the analysis symmetric case for which 
the components have equal masses of the particles m\ = ni2 = to, equal densities m — n 2 = no of 
the components and equal interaction constants 71 = 72 = 7- While this case is quite specific, it is 
possible to produce equal 7$ using a Feshbach resonance. The symmetric case may also correspond to 
two quasi two-dimensional Bose clouds with a strong dipole interaction separated by a rather high (but 
thin) barrier that suppress the tunneling (see 3 ) . On the qualitative level the results obtained for the 
symmetric case hold for the general case where such a symmetry between the components is broken. 
We assume the interaction between the particles of the same component is repulsive (7 > 0), and the 
stability condition with respect to a spatial separation of the components (7 > I712I) is fulfilled. Going 
ahead we note that in case of different velocities that condition is necessary but not sufficient one. 

If the temperature is much less than the temperature of the Bose-Einstein condensation the wave 
functions of the components can be presented as a sum of a large stationary part and a small fluctuating 
part 

^(r,t) = ^ 0i {r,t) + 6^i(r,t). (3) 
The stationary part of the condensate wave function reads as 

fe(M)=^e*( r 'e-¥, (4) 

in Ti 
where fa = — + (7 + 712)^0 are the chemical potentials of the components, and Sj = — Vy?j are 

their superfluid velocities. We will search for the fluctuating part of Eq. in the form: 

Sfl>i(T,t) = e*We"T kCrJe-*"* + vj(r)e*"**] . (5) 

Here the functions Uj and Vj are the plane waves 

Ui(r) = A ( e 4kr , v<(r) - B^ r . (6) 
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Substituting Eqs. © - j6]) into Eq. (J), in a linear in fluctuations approximation we obtain the following 
equation for the excitation energies hw and the eigenvectors: 



MV = uiV, 



(7) 



where 



( (k 



M 



1 + sik 



1 

T + l-sxk 



-7 



-7 

fc 2 

y + 1 + s 2 k 



-1 



2 



l-s 2 k 



(8) 



and 




Here and below all physical quantities are expressed in terms of dimensionless length and time 

Co 



t = 



-A 



(9) 



(10) 



where cq = y/ r yno/m is a sound velocity in a one-component condensate, and £ = fr/^/2mjno is a 
healing length. We also define the dimensionless parameter of the interspecie interaction 7' = 712/7. 

The u — v transformation procedure ([5]) is equivalent to the diagonalization of a quadratic form on 
Bose operators 18, 19|. The matrix M has four eigenvalues. The components of the eigenvectors of Eq. 
can be normalized as |Ai| 2 - |#i| 2 + \A 2 \ 2 - \B 2 \ 2 = ±1/V (V is the volume of the system). For 
two physical modes the norm of the eigenvector should be positive [19]. The spectra of the physical 
modes read as 



(si + s 2 )k 



+ 



\ 




(s„k) 2 + fc 4 7 ' 2 , 



(11) 



where s_ = Si — s 2 is the relative superfluid velocity. At si = s 2 = the spectrums (JTJ) has the 
standard Bogolyubov form 

uj± = k^/l± |y| +fc 2 /4. (12) 

As is clear from Eq. (|12[) . the stability condition with respect to a spatial separation 7 > I712I is the 
requirement for the excitation spectrum be real valued. A t small k the excitation spectrum is a sound 
one and the velocities of the modes are equal to c± = y/l ± The Landau criterium requires the 
spectrum (fTTj) be positive valued at all wave vectors. One can see from Eq. tfTTj) that at Si = s 2 = s 
the Landau criterium is reduced to the inequality s < c_. If only one component flows (s 2 = 0), the 
Landau criterium requires a fulfilment of the inequality s\ < c m = c_ •v/T + \"f\. If the velocity of a 
given component exceeds c m , the superfluidity condition is broken irrespective of a value and direction 
of the velocity of the other component [5], i.e. c m can be called the maximum critical velocity. 

At the velocities for which the energy ([TTj) is negative valued in some range of k the Landau criterium 
is broken and stationary waves can occur in the system. In contrast to a one-component system, in a 
two-component system an increase of superfluid velocities may result in that the spectrum of the lower 
mode be complex valued. As follows from Eq. (fTTj) . the spectrum remains real valued if the relative 
velocity satisfies the condition 

|si - s 2 | < c sop = 2c_. (13) 

At complex valued frequencies (fTTj) the amplitude of excitations grows with time that leads to a 
destruction of a homogeneous state and to a spatial separation of the components (or stratification). 



5 



In contrast to Ref . [7] , we consider the condition of stratification and the Landau criterium as different 
conditions. The Landau criterium yields a joint restriction on both superfluid velocities (and its mutual 
direction), while the condition of stratification is a restriction only on the modulus of the relative 
velocity. Note that c sop > c m , and under increase of the velocity of a given component, first, the 
Landau criterium is violated, and then, after further increase, the stability condition with respect to 
the stratification is broken. One can show, that the same situation takes place under simultaneous 
increase of two velocities (at nonzero relative velocity). Only at equal in value and oppositely directed 
superfluid velocities the Landau criterium and the stability condition are broken at the same point. 

Since we are interested in stationary waves in a homogeneous (not stratified) system, we will 
consider only the velocities for which the condition (fT3|) is satisfied. 



3 The stationary wave pattern 

Let us consider a two-component BEC that flows past an obstacle situated at the origin of coordinates. 
We assume the system is quasi-two-dimensional, i.e. it is thin enough to neglect the dependence of 
the condensate wave function on the transverse coordinate, and to consider all vector quantities as 
two-dimensional ones. If the size of an obstacle is much less than £ it can be considered as a point one. 
Under violation of the Landau criterium the obstacle behaves as a point source of waves (below we are 
only interested in stationary waves). The waves propagate from an obstacle with the group velocities 
v g ± = du)±/d\s., where cu± — ui±(k x , k y ) are given by Eq. (fTTj) at k z = 0. Here and below the index 
+ (— ) corresponds to the waves generated by the upper (lower) mode. The direction of the propagation 
for the wave with a given k is defined by the expression 

dbj±/dk v , 

tan * ± = ^k:' (14) 

where x+(x~) is the angle between the group velocity direction for a given mode and the axis x. 

For the stationary wave the frequency (fll]) is equal to zero, and the components of the wave vector 
k = {k x , k y ) are related by the equation 

Lu±(k Xl k y ) = 0. (15) 

It is convenient to use the angle 77 between the wave vector and the opposite direction of the axis x as 
an independent parameter, i.e. 

k = (— k COS77, k sin??). (16) 

Let us denote the angle between the velocities by 9 and select the axis x along the bisectrix of this 
angle (si = (s\ cos §, — s\ sin |) and s 2 = (s 2 cos |, s 2 sin |)). For the stationary waves Eq. ([15)) yields 
the following dependence of the wave number k on the angle 77 



-2 + si cos 2 (?? - |) + s 2 cos 2 (77+ |) T 

1/2 (17) 



\j^ a + (si cos 2 (77 — §) — s^cos 2 (t? + f))' 



A wave crest line is a line of a constant phase. For the stationary wave the phase can be obtained 
from the equation 



(r) = f kdr. (18) 
Jo 



In Eq. (|18p the integral is along a straight line going out from the origin of coordinates and directed 
parallel to the group velocity v g ± . The angle \i between k and r is defined by the expression 

/J = 7T-77-X± (19) 

(see Fig. H|). 
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According to Eq. (|18|) , the quantity r for the points at the wave crest with a given phase 4> satisfies 
the equation 

(20) 



k cos \x 

Substituting Eq. (fT9j) into Eq. (|20"|) . we get the equations for the stationary wave crest coordinates 
x = r cos x and y = r sin % in a parametric form 



x± = 



k± (77) cos 77 [1 — tan x± tan 77] ' 



= <?tanx± , 21 x 

k± (rj) cos 77 [1 — tan x± tan 77] 

In Eq. (|2ip the values tan%± arc the functions of the parameter 77. The explicit form of these functions 
is given by Eq. (|14[). in which after differentiation one should substitute Eqs. (fT6|) and (fT7l) . The range 
of values of 77 is determined by the condition for the function k±(rj) be real valued. Eqs. (|2ip allow to 
draw the stationary wave pattern for arbitrary values of s\, S2 and 9. 

Let us consider some special cases. For definiteness, we choose the parameter 7' = 0.5. Such a 
choice corresponds to c_ = 0.707, c m — 0.866 and c+ = 1.22 (in cq units). 

1. At equal superfluid velocities si = s 2 = s Eqs. fSPjl yield the expected result. At s > c_ a set 
of stationary waves corresponding to the lower mode arises. The crests for these waves are outside 
the cone bounded by the arms Xl — ± arcsin(c_/s). If the velocities s > c + are reached the second 
set of stationary waves appears. It corresponds to the upper mode and situated outside the arms 
X = ± arcsin(c+/s). Since s_ = 0, the stratification does not occur at any s. As an example, the 
stationary wave pattern for s = 1.5 is shown in Fig|2] 

2. If only one component (say, the component 1) flows, the stationary waves arise at the veloc- 
ities S\ > c m . These waves correspond to the lower mode. The waves are outside the cone X- = 
± arcsin(c m /s). In this case the frequency of the upper mode is always positive and the second set of 
waves cannot arise. At s% > c scp (for the parameters chosen c scp = 1.41) the system becomes unstable 
with respect to the stratification. The stationary wave pattern for the flow of one component with 
si = 1.0 is shown in FigO 

3. It is interesting to analyze the situation when the angle between the velocities is 9 = tt/2. In 
this case the obstacle emits waves only when the velocity of at least one component exceeds c m . Let us 
consider a more specific case of the velocities equal in magnitude s\ = S2 = s. Then the velocity range 
in which the stationary waves occur is limited by the condition c m < s < \[2c- (for the parameters 
chosen 0.866 < s < 1). In FigfJ]we present the stationary wave pattern for s = 0.9. One can see that 
in such a situation two sets of stationary waves arise. It is important to emphasize that both sets 
correspond to the lower mode (the frequency of the upper mode remains positive for the velocities in 
the range c m < s < \f2c-). 

In the limit 7' — > (that corresponds to the absence of the interaction between the components) 
each component should have its own set of stationary waves at any relative directions of the velocities. 
For rather large 7' this feature survives only in a close vicinity to 9 = tt/2. There should be a smooth 
transition from # = 7r/2to6' = 0. The analysis of the wave patterns at different 9 shows that at 
intermediate 9 a quite complicated pattern emerges: the crests of a given set end with cusps, and bridges 
connect them with crests from the other set. Under decrease of 9 the bridges and cusps disappear and 
the wave pattern becomes similar to one for 9 = 0. The stationary wave pattern with cusps and bridges 
is shown in Fig. [5] 

In general, the stationary wave pattern is qualitatively similar to one of presented in Figs. [2H 



4 The density pattern for the stationary waves 

There is a number of methods of probing BECs to get their density profiles (see [5D] ). Advanced technics 
was developed for the study of density profiles of spinor BECs 21,22,23,24,25] (that in certain sence 
can be considered as two-component ones). In particular, the density and the spin-density profiles of 
an atomic cloud were measured with a high resolution by the polarization-dependent phase-contrast 
imaging method [25U26j . 
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In view of modern experimental possibilities it is important to find the ratio of the components in 
the stationary waves and to determine specific features of the total density and the relative density 
patterns for the stationary wave in two-component systems. 

Since the stationary wave is the eigenmode, the ratio of the total density and relative density 
amplitudes can be found from the corresponding eigenvector (|9j). 

The density of a given component is the square modulus of the condensate wave function pi = ■ 
This quantity can be presented as a sum of the unperturbed density no and the perturbation Srii caused 
by the stationary wave. Using the equation for the condensate wave function ([3]) and Eqs. (J4j) — dSJ) , 
we obtain 

Srii = ip Sip* + S^i/jq = 5p % cos(kr), (22) 

where 



5p i = A^Z{A i + B i ). (23) 

here A is the amplitude of a given eigenmode. The amplitude A depends on the intensity of the source 
(obstacle) and on the distance from the source. As follows from Eqs. (|7I8[) , the components of the 
eigenvectors are real valued quantities. That is why Eq. (|23j) may correspond to in-phase oscillations 
of the densities or to the oscillation with the phase shift equal to ir (out-of-phase). Here we consider 
the case 7' > for which the lower Bogolyubov mode corresponds to out-of-phase oscillations and the 
higher mode - to in-phase oscillations. 

The total density and the relative density oscillation amplitudes in the stationary wave are Sp± = 
pi ± pi. The specific of the symmetric case (components with equal masses, densities and interaction 
constants) is that in the absence of the flow the oscillations of the total density vanish for the lower 
Bogolyubov mode (8p+ = 0), and the oscillations of the relative density are absent in the higher mode 
(5p- = 0). As was shown in the previous section, in most cases the stationary waves are caused by the 
lower mode. Therefore it is important to clarify whether is the total density disturbed in the stationary 
waves. 

Using the analytical expressions for the eigenvectors of the matrix |8]) one finds that the exact 
relations Spi ~ ~Sp2 for the lower mode and Spi — 8p% for the upper mode are hold in case of equal 
velocities (si = S2). In means that in the stationary wave pattern shown in Fig. [2] the lower mode set 
is visible in the relative density image, while the upper mode set - in the total density image. 

In the general case S\ ^ S2 we find that Spi ^ —5p2 and the total density oscillations are nonzero 
for the stationary waves that correspond to the lower mode, in particular, for the pattern shown in 
Figs.[3]-[5l For these cases the ratios | <5/^i | / 1 (5/02 1 along the crests are shown in Fig. [6] Since this ratio 
differs from unity, both, 5p+ and 5p- are nonzero and stationary wave patterns should be visible in 
the total density image as well as in the relative density image. We also would like to point out the 
feature that follows from Fig. [6] If only one component flows past the obstacle the stationary waves 
contain admixture of both components, but the flowing component has larger amplitude. 

In cases shown in Figs. 13141 two or three sets of stationary waves interfere. Thus the density pattern 
depends not only on the eigenvectors but on the amplitudes of the eigenmodes, and another approach 
should be used to analyze the density profiles. Here we use the approach developed in [17] . 

The interaction of the Bose gas with the obstacle is described by the Hamiltonian 



H in t = J dr^ViO-Jhfc 



where Vj (r) is the potential of the interaction between the i-th component and the obstacle. Considering 
the obstacle as a point source and assuming that the obstacle interacts identically with both components 
we set Vi(r) = V^r) = VoS(r). Respectively, the terms Vo5(r)ipi should be added to the right hand 
parts of the Gross-Pitaevskii equations |T]). The Gross-Pitaevskii equations can be rewritten in terms 
of density and velocity fields (the hydrodynamic form). Linearizing the hydrodynamic equations with 
respect to the density and velocity fluctuations (see, for instance, [5]), and excluding the velocity 
fluctuations we arrive at the equations for the Fourier components of Srif. 



is t k) 2 + £ ; 2 (l + ^ 



fei(k) + fcV5n 3 -i(k) = -k 2 V , (24) 
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where Vq = Vq/j. In deriving Eq. ([24]) we take into account that all time derivatives are zero for the 
stationary waves. 

Solving Eq. ([24]) and taking the inverse Fourier transformation we get 
Ay /■ 7r / 2 r°° \ _ <yi i hL _ ( s i k ) 2 

]_„,,/" I * vV(*-~g) (25) 

where the quantities fc± are given by Eq. (fT7)l . The angle is determined by Eqs. I|19[) and (HH). 
The integral over fc in Eq. (|25| can be evaluated analytically with the use of the residue theorem. The 
poles k — and k = k + (for real k± ) yield the contribution of the lower and the upper Bogolyubov 
modes, respectively. The integral over rj is evaluated by the stationary phase method. The number of 
stationary phase points coincides with the number interfered waves in the stationary wave pattern (1 
point for Figs. [2] and [3l 2 points for Fig. [4] and 3 points for Fig. [5]). 
The answer can be presented in the following form 

Ci,x X 7= > (26) 

where the index A numbers the sets that contribute to the stationary wave pattern in a given sector 
of x- I n Eq. (f2"rJ|) short range terms oc 1/r 2 are omitted. The coefficients C^xfc) an d n i,\{x) are quite 
complicate and we do not present the explicit expressions here. 

The case si — S2 = s can be analyzed directly from Eq. ([23)1 . At such velocities 



5n ' = 6n " = ^e^j_^j Q kdk w^> (27) 



where k + = 2y s 2 cos 2 -q — c 2 + . One can see from Eq. ([27]) that, first, 8n\ — 5ri2 and the relative density 

remains unperturbed, and, second, the stationary waves patterns are caused by the upper mode only 
and emerges at s > c+. Such features have clear explanation: the obstacle with V\ — V2 cannot excite 
the mode with 5ri\ = —5ri2. We emphasize that these features are specific for the symmetric two- 
component condensate. If the bare modes of the components differ from each other the lower mode is 
not pure relative density oscillations (Sni =/= —Sns) and the obstacle with V\ = V% may excite the lower 
as well as the upper mode. 

Using Eq. (|26[) and explicit expressions for Ci y \ and k^a we obtain the density plots for the interfered 
stationary waves given in Figs. [H The results for the total density and relative density patterns are 
shown in Figs. [3 [H One can see that the images presented have clear interference structure and both 
the total density and relative density measurements can be used for the visualization of the stationary 
waves. Note that in the case considered the relative density pattern in much more contrast than the 
total density pattern. 

Ending this section we would note the following. The excitation of stationary and propagating 
waves by the obstacle can be considered as a kind of Cherenkov radiation. Cherenkov radiation arises 
when the velocity of a radiating object in the medium exceeds the phase velocity of radiated waves in 
this medium. In this respect it was unclear why the critical velocity may be larger than the velocity of 
the lower mode c_. The answer is the following. At relative motion of the components the eigenvector 
for the lower mode differs from one for the condensate at rest or at s_ = 0. In other words, in the flow 
with s_ 7^ the structure of the modes is modified, and the mode radiated at nonzero relative velocity 
is not the c_ lower Bogolyubov mode. 



5 Conclusions 

In conclusion, the properties of stationary waves arising in a flow of a two-component BEC past an 
obstacle have been studied. The problem was considered for a special case of the symmetric two- 
component system (with the components of equal masses of the particles, equal densities and equal 
interaction constants). Nevertheless, the majority of conclusions are applicable also to the cases in 
which such a symmetry between the components is absent. Let us recite these conclusions. 
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Fig. 1 The scheme of a wavefront crest 

In a two-component flowing superfluid system the energies of the excitations can take on not only 
negative, but complex values. At reaching of negative values the Landau criterium is broken, and at 
reaching of complex values the system becomes unstable with respect to a spatial separation. The 
Landau criterium is a joint condition of both superfluid velocities (and the angle between them) in the 
lab reference frame. The stability condition with respect to a spatial separation is a condition solely 
on an absolute value of the relative velocity (which does not depend on the reference frame). In case 
of equal (in modulus and in direction) velocities of the components a spatial separation does not arise 
at any velocities. Under increase of the velocities the Landau criterion is broken first, and then the 
stability condition is broken. Stationary waves arise when the Landau criterion is already broken, but 
the system remains stable with respect to a spatial separation. At equal velocities the stationary waves 
generated by the lower and the upper modes can arise. If only one component flows, the one set of 
stationary waves (corresponding to the lower mode) may emerge. If the angle between the velocities is 
close to 7r/2, two sets of stationary waves can arise, and both of them correspond to the lower mode. 
In general, the stationary waves are visible at the total density and relative density images. 

We did not consider here the ways of creation of relative flow of the components in two-component 
atomic vapors, but we would like to mention some other possibilities. It is quite simple to realize such 
a flow for two components separated with a thin barrier. In this case the interaction between the 
components should contain a long-range (for example, dipole) part. Strictly speaking, in such systems 
the spectrum of excitations may differ from (jlip (due to a long-range interaction). Nevertheless one 
can expect that the stationary waves will behave qualitatively the same as in the case considered in 
this paper. Similar phenomena may emerge in some other systems, for example, multilayer electronic 
systems with superfluid indirect excitons [27]. Another perspective object is superfluid polaritons in 
semiconductor microcavities, where superflow can be controlled by the laser beam. In a one-component 
polariton system the observation of stationary waves was reported recently [28] . Relative motion of the 
components may also arise if one component is electrically charged and the system is subjected by an 
electromagnetic field. For instance, such a situation takes place in neutron stars [29 . 
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Fig. 2 The wave pattern forming at the flow of components with equal velocities (s = 1.5). The waves caused 
by the lower and the upper modes are shown with solid and dashed lines, respectively. 




Fig. 3 The wave pattern forming at the flow of one component with the velocity Si = 1.0 (s2 = 0). The axis 
x is along the flow direction. 




Fig. 4 The wave pattern forming by two orthogonally directed flows with the velocities s 1 — 
directions of the velocities are shown by the arrows. 
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Fig. 5 The wave pattern at si = S2 = 1 and # = 0.45-7T 
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Fig. 6 The ratio of densities of the components in the stationary waves along the chests for: a) one component 
flow (si = 1, S2 = 0); b) two component flow in mutually orthogonal directions (9 — tv/2, si = S2 = 0.9); c) the 
case of 6 = 0.45-7T, si = S2 = 1. 
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Fig. 7 The total density (a) and the relative density (b) patterns (in relative units) for the stationary waves 
at 6 = 7r/2 and si — S2 = 0.9. 
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